Non-Gaussianity and large-scale structure in a two-field inflationary model 
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Single field inflationary models predict nearly Gaussian initial conditions and hence a detection of 
non-Gaussianity would be a signature of the more complex inflationary scenarios. In this paper we 
study the effect on the cosmic microwave background and on large scale structure from primordial 
non-Gaussianity in a two-field inflationary model in which both the inflaton and curvaton contribute 
to the density perturbations. We show that in addition to the previously described enhancement 
of the galaxy bias on large scales, this setup results in large-scale stochasticity. We provide joint 
constraints on the local non-Gaussianity parameter /nl and the ratio £ of the amplitude of primordial 
perturbations due to the inflaton and curvaton using WMAP and SDSS data. 
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I. INTRODUCTION 

One of the most important questions that cosmology 
faces today is the origin of structure in the universe. 
The generally accepted paradigm is that of inflation [J-Q 
which produces small adiabatic perturbations that evolve 
into the observed structure. The inflationary paradigm is 
extremely powerful as it remedies most of the problems 
of the original Big Bang scenario and also has a set of 
predictions that are well confirmed by current observa- 
tions. On the other hand, although the generic predic- 
tions of inflation are quite clear, the nature of specific 
physical processes that govern inflation are still poorly 
understood. 

The major obstacle in understanding inflation is that 
it can not be directly observed either in the laboratory 
or with telescopes. This problem is at the same time a 
virtue of inflation as it allows to indirectly probe physics 
at energies and time-scales that are far beyond the reach 
of current facilities. By comparing astrophysical observa- 
tions with predictions of various inflationary models one 
can expect to distinguish between different extensions of 
the Standard Model of particle physics ||. Understand- 
ing of the reheating phase of inflation can provide a link 
between scalar fields driving inflation and the observable 
Universe that consists of dark and baryonic matter. 

One of the many possible ways to deeper under- 
stand inflation is by studying the primordial density 
fluctuations. The usual inflationary model of a slowly- 
rolling inflaton field requires that the perturbations are 
highly Gaussian (fHH and hence the detection of non- 
Gaussianity in either the cosmic microwave background 
(CMB) spectrum or the large scale structure (LSS) distri- 
bution would be a clear evidence that the physics driving 
inflation is more complicated than the standard inflaton 
scenario. 

Non-Gaussianity naturally arises in inflationary mod- 
els with more than one field |10l4l3| . One of the most 
studied models is the curvaton model [l3l - [l8j ]. in which 
initial perturbations are generated by the curvaton field 



after inflation is over. In this model significant non- 
Gaussianity can be generated since the predicted curva- 
ture perturbation is proportional to the square of the cur- 
vaton field (as distinct from single-field inflation, where 
the required smoothness of the inflaton potential renders 
the curvature perturbation very nearly linear in the field 
fluctuations). 

Most attempts to constrain non-Gaussianity have used 
the so-called "local-type" or /nl parameterization [l9| in 
which one includes a quadratic term into the primordial 
potential, $ = 4> + /nl^ 2 - In this parametrization both 
linear and quadratic terms in the potential originate from 
the same Gaussian field, e.g. a curvaton field, and the 
contributions from perturbations in other fields (e.g. the 
inflaton field responsible for inflation itself) are negligi- 
ble. 

The signature of local-type non-Gaussianity in the 
CMB has been described at length [20(. It has also been 
established that /nl has an effect on the galaxy bispec- 
trum [2l| - [23j . The effect on the large-scale gala xy p ower 
spectrum has been considered only recently [2J-[28|, but 
it rapidly became clear that the method was competitive, 
stimulating work on iV-body simulations of halo forma- 
tion in non-Gaussian cosmologies [29l432| . Recent con- 
straints have been derived from the CMB bispectrum as 
measured by WMAP (33T - [39f and from large scale struc- 
ture in the Sloan Digital Sky Survey (SDSS) [H|. Re- 
cently, ~ 3cr evidence for excess clustering consistent with 
non-Gaussianity has been identified in the NRAO VLA 
Sky Survey (NVSS) 0. 

In this paper we extend the formalism to include the 
case where both the inflaton and curvaton contribute sig- 
nificantly to the curvature perturbation. Perturbations 
generated by the inflaton field are purely Gaussian, while 
curvaton fluctuations can result in non-Gaussianity if the 
conversion from curvaton fluctuation 8a to primordial 
potential $ contains quadratic terms. The ratio of infla- 
ton to curvaton contributions £ is arbitrary: the frame- 
work of the curvaton model allows it to take on any posi- 
tive value. Usually one takes £ ^> 1 since in the opposite 
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limit <C 1) the curvaton has no observable effect on the 
primordial perturbations. Here we investigate the conse- 
quences of general £ - including values of order unity - for 
the CMB and LSS. The type of non-Gaussianity gener- 
ated could be called "local-stochastic," in that it results 
from local nonlinear evolution of the inflaton and curva- 
ton fields (and thus the primordial bispectrum will have 
the local- type configuration dependence), but that the 
full nonlinear potential $ is not a deterministic function 
of the linear potential. 

Studying non-Gaussianity is particularly important in 
the face of the current generation of CMB projects 
such as the Planck satellite as well as ongoing and future 
LSS projects. To fully exploit the potential of the future 
probes it is imperative to investigate theoretically the 
range of types of non-Gaussianities that can be produced 
in unconventional inflation (e.g. multi-field models), and 
understand what effect they have on the CMB and LSS. 

The rest of the paper is organized as follows. In Sec. ITU 
we discuss the generation of non-Gaussian primordial 
perturbations in the inflationary model with both infla- 
ton and curvaton fields contributing to the curvature per- 
turbation. In Sec. Mil we describe the effect of two-field 
models on the CMB bispectrum. In Sec. II VI we derive the 
halo power spectrum using the peak-background split for- 
malism (4^ . and in Sec. [V] we consider the angular power 
spectrum of galaxies. Section PVTl provides the constraints 
on the two-field model from existing data, and we con- 
clude in Sec. EE 



II. NON-GAUSSIAN INITIAL 
PERTURBATIONS IN TWO-FIELD 
INFLATIONARY MODELS 

We consider a model of inflation where both the infla- 
ton and the curvaton contribute to the primordial density 
perturbations. This configuration can exhibit a rich set 
of phenomenology, including both non-Gaussianity and 
various mixtures of adiabatic and isocurvature perturba- 
tions [43M47| . In this paper, we will restrict ourselves 
to the case where the dark matter decouples after the 
curvaton decays and its energy density is thermalized. 
This ensures that no dark matter isocurvature perturba- 
tion is produced, and the only observable perturbation is 
the curvature perturbation £ that is conserved between 
curvaton decay and horizon entry. 

The simplest case is that of two non-interacting scalar 
fields: the inflaton ip and the curvaton a. The latter is 
taken to have a quadratic potential, 



1 



= - m a . 



(1) 



During inflation, the inflaton dominates the energy den- 
sity of the Universe, whereas the curvaton is effectively 
massless (m <C H) and pinned by Hubble friction to a 
fixed value a (aside from perturbations to be described 
later) . Quantum fluctuations generate a spectrum of per- 



turbations dip and da in both the inflaton and curvaton 
fields: 
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and Psa{k) = 



2k 3 ' 



(2) 



where i?* is the Hubble rate evaluated at the horizon 
crossing for a given mode, i.e. when k — aH, and the 
post-horizon-exit field perturbations arc defined on the 
uniform total density slice. The if and a perturbations 
are nearly Gaussian and uncorrelated. 

The inflaton perturbation is parallel to the unper- 
turbed trajectory in (tp, <r)-space and hence is an adia- 
batic perturbation; indeed it behaves the same way that 
perturbations behave in single-field inflation. The cur- 
vaton perturbation however is an isocurvature perturba- 
tion and can have complicated dynamics. In the simplest 
version of the curvaton scenario, the curvaton begins to 
oscillate after the end of inflation when the Hubble rate 
drops to H ~ m. As a massive scalar with zero spatial 
momentum, its energy density subsequently redshifts as 
p a oc a~ 3 . Of interest to us is the fact that for quadratic 
potentials (Eq. [TJ) this energy density is also proportional 
to the square of the curvaton field a = a + Sa, i.e. 



5 _ <w = 

Pa a 



Sa 2 - {5a 2 ) 



(3) 



where the subtraction of the variance arises from the 
0(Sa 2 ) expansion of p a oc a 2 + (5a 2 ). The quadratic term 
allows the curvaton to generate a non-Gaussian density 
perturbation. In the radiation-dominated era, the cur- 
vaton's contribution to the energy density increases as 
Oct oc a, thereby enhancing the importance of 8 a . The 
decay of the curvaton and the thermalization of its energy 
density result in a non-Gaussian adiabatic perturbation. 

The SN formalism [48[ extended into the nonlinear 
regime Q quantitatively provides the curvature pertur- 
bation to second order in the field perturbations (dtp, Sa); 
this is (Eq. 26 of Langlois et ai. [461 ] ^ 
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where the subscript » denotes evaluation at horizon exit, 
and r is related to the fraction of the energy density in 
the curvaton when it decays: 
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4/9rad + 3/9ct 



(5) 
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The primordial potential perturbation in the Newto- 
nian gauge is then given by the usual expression $ = 
|£ (e.g. [49]; but note that in large scale structure 
non-Gaussianity studies the opposite sign convention is 
adopted, so that $ > corresponds to overdensities). 

We may put the primordial potential in a form more 
closely related to that of large-scale structure non- 
Gaussianity studies: 



$(x)=^(x)+0 2 (x) + / NL [^(x) 



(G) 



3 



where <j>\ and </> 2 are the parts of the linear primordial po- 
tential corresponding to the inflaton and curvaton fields 
respectively. Their power spectra are given by 



k 3 9 / Hi 



and 



^(*) 



4r 2 



25 \2tt(t, 



The non-Gaussianity parameter is 



(7) 



(8) 



(9) 



(We use the tilde since the label "/nl" is usually used to 
denote the non-Gaussianity parameter appearing in the 
primordial bispectrum.) 

It is convenient to specify the relative contribution 
of the inflaton and curvaton fields to the primordial 
potential using the ratio of standard deviations £ = 
cr(4>\)/a(4>2)- Thus a fraction £ 2 /(l + £ 2 ) of the power 
comes from the inflaton and a fraction 1/(1 + £ 2 ) from 
the curvaton. This ratio is 
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{d(f>/dN)* 



(10) 



where N is the number of e-folds remaining in inflation. 
Thus the observable features of this model are specified 
by the primordial power spectrum P$(fc) and by the two 
new parameters /nl and £ (in principle £ will have a scale 
dependence d In £/d In k of order the slow roll parameters, 
but unless non-Gaussianity is detected at high statistical 
significance this cannot be measured). We will work with 
these parameters from here forward. 



III. THE CMB BISPECTRUM 

The effect of local-type non-Gaussianity on the CMB 
bispectrum has a long history, both in purely adiabatic 
models as considered here, and in locally non-Gaussian 
isocurvature models (HO, HU ■ We evaluate the CMB bis- 
pectrum using our set of parameters here. 

The CMB constraints on primordial non-Gaussianity 
come from the measurements of the CMB angular bis- 
pectrum fTgj j . 



(11) 



where a£ rn is the CMB temperature fluctuation expanded 
in spherical harmonics: 



a lm ee / d 2 n^By/ m (n). 



(12) 



If the primordial fluctuations are adiabatic scalar fluctu- 
ations, then ai m can be easily expressed in terms of the 



primordial potential $ and the radiation transfer func- 
tion ge(k): 



aim = 47r(-i) 



d 3 k 

(2tt) 3 



$(k) 5 ,(fc)Y/ m (k). (13) 



From Eq. ^ it follows that in the Fourier space primor- 
dial potential can be decomposed into parts associated 
with the linear potential perturbations <j>\ and (fa, and 
with the nonlinear coupling /nl: 

$(k) =0i(k)+0 2 (k)+0 WL (k). (14) 

Here </>jvx(k) is the /NL-dependent part, 

d 3 p 



4>nl0<) = /nl 



(2tt) 3 



,(k + p)0*(p). (15) 



Using Wick's theorem we calculate the bispectrum of the 
total potential, P$(fci, fc 2 , ^3). It contains one contribu- 
tion from allowing each of the $(kj) to have a contribu- 
tion from /nl; in the case where this is £3: 

([01 (ki) + 2 (ki)][<Mk 2 ) + 2 (k 2 )]0Ar L (k 3 )) = 

2(2ir) 3 5^(k 1 + k 2 + k 3 ) - hh P^P^ki), (16) 

(1 + 4 ) 

where P$(fc) is the power spectrum of the total potential. 
The total bispectrum is the sum of this and the similar 
contributions where /nl contributes to hi and to k 2 : 



/ 



B*(fci,fa,*3) = 2-^WP$(fc 1 )P$(fc 2 ) + 
(1 +CT 

PMPzih) + P*(fca)P*(*3)]. ( 17 ) 

For constraining non-Gaussianity it is convenient to 
introduce two new variables: x\ = /nl/(1 + £ 2 ) 2 and 
X2 = 1/(1 + £ 2 ). Using the bispectrum of ^(k), we can 
finally write the CMB angular bispectrum (via a calcu- 
lation similar to Ref. [Hj|): 



nmim 2 m3 
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r dr 



xK(r)bi(r)bZ h (r) + l£(r)C {r)b^{r) 



where G™\™£ m3 is the Gaunt integral, and we use 



bf(r) = - k 2 dkP^(k)g(k)j e (kr) (19) 
Jo 



and 



&T(r) 



2xi 

7T 



k 2 dkg{k)ji{kr) 



(20) 



We see from Eq. (fl~7|) that CMB bispectrum measure- 
ments of /nl that assume a purely curvaton contribution 
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and short- wave length pieces: 




FIG. 1: The allowed range (2a) of /nl as a function of X2 
derived from the WMAP data [3q. As discussed in the text, 
/nl becomes unconstrained as X2 — ► because in this case 
the statistics describing the density distribution is dominated 
by the inflaton field and is nearly Gaussian. 



p(x) =p(l + Si+5 s ). 



(21) 



The linear density field here is a sum of two independent 
Gaussian components S = #W _|_ j( 2 ) originating from the 
inflaton and curvaton fields, respectively. 

The local Lagrangian density of halos now depends 
on the large scale matter perturbations of both gaussian 
fields: 



n(x) = n[l + b\5 : 



(i) 



?( 2 )l 



(22) 



We will see that bi can be scale-dependent (fc-dependent), 
which means that in position space it should be thought 
of as a (possibly nonlocal) operator acting on the density 
field, i.e 

Fourier transform of 62 (k) 

We can express the bias parameters in terms of the 
number density function 



62 <5 Z (x) is the convolution of 5) (x) with the 



On 



OS 



(23) 



(£ = 0) are actually measuring xi in the more general 
case. 

Clearly, when the curvaton field dominates the pertur- 
bation power (x2 ~ 1), we have /nl = /nl, and our 
/nl constraints are identical to the constraints on mod- 
els with negligible inflaton perturbation. However, as the 
contribution from the inflaton field increases (#2 — > 0) the 
role of the contribution of the curvaton field to the pri- 
mordial curvature perturbation becomes negligible and 
the statistics describing the density distribution becomes 
very nearly Gaussian. In that case /nl becomes com- 
pletely unconstrained (Fig. Q]) and we can no longer make 
robust predictions regarding the presence of the second 
field on the basis of the non-Gaussianity observations 
alone. The CMB bispectrum is therefore not capable 
of breaking the degeneracy between /nl and £. The use 
of other constraints is necessary. In this paper, we will 
use large-scale structure, although we note that in prin- 
ciple the CMB trispcctrum might also be useful for this 
purpose, since it scales as / NL /(1 + £ 2 ) 3 = x\/%2 and 
thus in combination with the bispectrum would allow the 
(/nl,£) degeneracy to be broken. 



It is easy to check that the bias for the 8\ field is just 
the usual Lagrangian bias that applies to Gaussian cos- 
mologies; for example, in the Press- Schechter model [52[ 
it is b\ = b g = 5 c /afi — S^ 1 with S c = 1.686 quantifying 
the spherical collapse linear over-density. To calculate 62 
we note that short-wavelength modes S s in an overdense 
region determined by 81 can be written as: 



(1 + 2.wPV1 2) +/nl(^ 2) ) 2 + ^ 1) " 

= a\x 1< p^+X 2 (^) 2 + ^] , (24) 



where X x = 1 + 2/ NL 0, (2) and X- 



NL- 



Here a is 



the transfer function that converts the potential into the 
density field, 5(k) — a(fc)$(fc). In general one may think 
of a as an operator defined by its action in Fourier space, 
i.e. when applied to a real-space function such as </>(x), 
we have 



<fye- ik '4,(y). (25) 



The specific function a(k) is given by Eq. (7) of Slosar 
et al. 1251: 



IV. LOCAL NON-GAUSSIANITY IN 
PEAK-BACKGROUND SPLIT FORMALISM 



In this section we outline the inflationary scenario in 
which both inflaton and curvaton fields are contribut- 
ing to the initial curvature perturbation and derive the 
power spectrum of galaxies using the peak-background 
split formalism (42[. 

We decompose the density field into a long-wavelength 



a(k; z) 



2c 2 k 2 T{k)D(z) 



30, 



H 2 



(26) 



where T(k) is the linear transfer function with conven- 
tional normalization T(0) = 1, and D{a) is the growth 
function normalized to D(z) = (l + z) -1 at high redshift. 
The inverse operator a -1 is obtained by the replacement 
a(k) -> a~ l (k) = l/a(k). 

This shows that local number density in the non- 
Gaussian case depends not only on Si, but also on X\, 
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Xn, and hence b 2 becomes 



&2 = 



On 



_c^ (2) (x) 



_ 1 dn 



dXi 



(27) 



where the derivative is taken at the mean value X\ = 1. 

We can further simplify this expression by considering 
a rescaling of the power spectrum on the small scales due 
to the presence of non-Gaussianity. In a "local" region of 
some size R, and for small-scale Fourier modes k ^> R^ 1 
within this region, there is a local power spectrum 



P. 



local 



(k) = 



e + (l + 2f NL a-^ 



X Pf oh *\k), 



from which we obtain the rescaling of as : 



_local 



(28) 



(29) 



Using these expressions we can change the derivatives in 
Equation (|27|) to finally obtain 



b 2 (k) 



2/i 



NL 



a-^k) 



d Inn 
dlno-*° ca1 ' 



(30) 



For further calculations we assume the mass function 
to be universal, i.e. we assume that it depends only on 
the significance function v(M) = 8 2 /a 2 (M): 



n{M, v) = M- 2 vf(v) 



dbxM' 



(31) 



This assumption is much more general than the Press- 
Schechter picture, e.g. it holds for the Sheth-Tormen 
mass function (53l | as well. Universality implies that 



d Inn 
d In cr' ocal 

from which we derive 



S c b g , 



b 2 (k)=b g + ^^-b g a- 1 (k). 



(32) 



(33) 



The standard Gaussian bias in Eulerian space is given 



by b 



1. The halo overdensity in Eulierian space 



in the non-Gaussian case is then obtained by multiplying 
Eq. by 1 + (5 ; ; to first order, 



n(x) 



Si + + b 2 5f ] 



[l+hmsw 



[l + b 2 (k)]Sl 2) . 



(34) 



We can now write down the halo power spectrum in the 
form: 



Phh(fc) 



(l + 6i) 2 g 2 + [l + & 2 (fc)]' 



P hn (k). (35) 



Finally, plugging in bi(fc) 
b 2 (k), we obtain 



b and using Eq. (f33|) for 



Phh(fc) 



eb 2 + [b+ 2(6 - i)/niA(i + er^-^k)] 2 



i+e 



xP lin (fc). 



(36) 



In the limit of ^ — > 0, i.e. when the contribution from the 
inflaton field is negligible we recover the standard expres- 
sion for the power spectrum with the curvaton generated 
non-Gaussianity (Eq. 32 of (25j). 

(It should be noted that the non-Gaussianity also in- 
troduces small corrections to the scale-independent part 
of the bias, because the small-scale fluctuations that must 
collapse to form a massive halo have a non-Gaussian den- 
sity distribution. This effect has been seen in simula- 
tions with /NL-type non-Gaussianity, where it is negative 
for /nl > 0, resulting in a slight reduction of the non- 
Gaussian bias enhancement at large k, and even a change 
in sign of the /nl effect at very small scales 0,|3l|. How- 
ever, since current studies of non-Gaussianity using LSS 
allow the scale-independent bias to be a free parameter, 
they are not sensitive to this effect; it would only be 
important if the Gaussian contribution to the bias were 
inferred indepedently, e.g. from measurements of halo 
mass and the mass-bias relation.) 

A further consequence of this model that does not arise 
in the case with only curvaton perturbations is large- 
scale stochasticity. In particular, the squared correlation 
coefficient 



X(k) 



^hm(^) 



Phh(k)P mm (k) 



(37) 



deviates from unity on the largest scales. We can see this 
by writing the cross power spectrum Phm(fc) as 

pUk) = ii±M^il±MMp lin(fc) . (38) 

1 + £ 

In the linear Gaussian theory, one would have x = 1) 
whereas in our case we have 



{(i+b^e + ii+hjk)}} 2 
(i+e){(i+bi) 2 e+[i+Hk)?y 



(39) 



Note that if /nl ^ 0, on large scales |&2(&)| ^ 1, t>i and 
hence 



(40) 



An example of the onset of scale-dependent bias and 
stochasticity is shown in Figure [5] note that this type 
of stochasticity effect exists only for x 2 1. 

It is important to note that stochasticity can arise even 
for x 2 — 1 in two ways. One is that on small scales, there 
is a breakdown of linear biasing. However, since our con- 
straints on non-Gaussianity arise from the largest scales 
in the survey (mainly the I < 25 quasar data points) this 



6 



effect can be neglected. The other is halo shot noise (e.g. 
54]), which arises from the fact that haloes containing 
multiple galaxies (or quasars) can produce a large "1- 
halo" contribution to the correlation function at small 
separations. When transformed to Fourier space at large 
scales (small k), this results in additional contribution to 
the power spectrum of 



Bias 



lim P 1 haio(fc) 



/ 



47rr 2 £i halo (r) dr, (41) 



where £i hak>(?") is one- halo correlation function. In prin- 
ciple since P(k) oc k on large scales, the halo shot noise 
term (oc A; ) may become important; since it is random 
and not determined by the underlying long-wavelength 
modes of the density field, it also produces stochasticity. 

However, the halo shot noise is expected to be a very 
small contribution for our data sets. A simple way to 
see this is to note that the ratio of the halo shot noise 
to the usual shot noise is equal to twice the ratio of 1- 
halo pairs to the number of galaxies (this follows from 
Eq. (|4"Tj) and the definition of the correlation function). 
For the quasar sample, a simple counts-in-pixels analysis 
of the catalog suggests that 0.6% of the quasars are in 
pairs (the Healpix pixels |55j used are 3.5 arcmin in size, 
i.e. much larger than haloes at z > 1), suggesting that 
the halo shot noise term is Cj.ihaio ~ 0.012/tt,2d- This is 
two orders of magnitude smaller than the error bars on 
the lowest-/ quasar autopower point displayed in Slosar 
et al. [25[ and hence negligible. 



V. GALAXY ANGULAR POWER SPECTRUM 

Additional constraints on the primordial non- 
Gaussianity come from the observations of the Large 
Scale Structure (LSS), and in particular from large 
galaxy surveys such as SDSS. These constraints are pri- 
marily driven by extremely large scales, corresponding to 
wave numbers k < O.OlAfpc" 1 . The data sets used for 
LSS constraints include spectroscopic and photometric 
Luminous Red Galaxies (LRGs) from SDSS, photometric 
quasars from SDSS and cross-correlation between galax- 
ies and dark matter via Integrated Sach- Wolfe effect. De- 
tailed descri ptio n of the data used for LSS constraints can 
be found in |25l . |56| and here we will only emphasize the 
redshift ranges of the most important datasets we used. 

The photometric LRGs dataset was constructed and 
discussed in detail in jszj. The sample was sliced into 
two redshift bins with 0.2 < z < 0.4 and 0.4 < z < 0.6. 
Our spectroscopic LRG power spectrum comes from (58j 
and is based on a galaxy sample that covers 4000 square 
degrees of sky over the redshift range 0.16 < z < 0.47. 
Quasars for our constraints were photometrically selected 
from the SDSS quasar catalog consisting of UVX ob- 
jects and DR3 catalog !6Qj. Quasars in our sam- 
ple fall into two redshift bins with 0.65 < z < 1.45 and 
1.45 < z < 2.0. 
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FIG. 2: The bias and stochasticity of galaxies at z = 1 in a 
model with x\ = 30 and X2 = 0.5 (/nl = 120). The solid 
lines show a tracer with 6 = 2 and the dashed lines a tracer 
with b = 3. The background cosmology and power spectrum 
are those of WMAP5. 



The tightest LSS constraints on primordial non- 
Gaussianity involve purely photometric samples (where 
one observes the two-dimensional projection of the galaxy 
distribution with a set of color cuts) since this allows the 
largest volume to be probed with the highest number 
density. At low k, where the effects of primordial non- 
Gaussianity on the power spectrum are largest, there is 
less of an advantage to having the large number of modes 
(oc &max instead of achievable via spectroscopy. In 

order to obtain constraints on the /nl and £ parameters 
from LSS we need to modify Eq. (|3l)|) to give the angular 
power spectrum. 

To obtain the angular power spectrum, we project the 
galaxy density field 8 g ^D along the line of sight n and 
take into account redshift distortions I61rl63|: 



l + 5 g (n) = J dy f(s)[l + 5 g , 3D (y,yn) 



(42) 



7 



Here, s = y + H~ 1 v ■ n is the redshift distance, f(s) 
is the normalized radial selection function, and we have 
explicitly written the mean contribution to the density 
field. We further note that peculiar velocities are gener- 
ally small compared to the size of the redshift slice and 
hence we can Taylor expand selection function as: 

f(s)^f(y) + ^-^-v(yn)-h. (43) 
all ay 

At large scales where the density perturbations are <C 1, 
we may ignore second-order terms, i.e. we may ignore 
the product of the velocity term in Eq. (|4"3"]) with the 
S g ^i}(y,yh) term in Eq. (|4"2"j) . This allows us to split S g 
into two terms as S g = 6® + S g . In terms of the Fourier 
transformed fields, we can write S g and 8 r g as: 



and 
S r g (n) = 



dyf(y) J -^^8 g , 3D (y,k)e- 
df f d 3 k 1 



ik-ny 



(44) 



d v~r 



dy J (2tt) 3 aH 



v(y,k) -ne-*' 6 " . (45) 



The velocity can be related to the density perturbation 
using linearized continuity equation: 



H- x v(y,V 



0.6 



k 

'jfe 5 



^S g (y,k)-. (46) 



We can further transfer redshift dependence of 6 g into 
a growth function D(y) and expand S's in Legendre poly- 
nomials Pe(x) using the following identity: 



e - lk .ny =J2(2t + l)i e j e (ky)P e (k • n) 



(47) 



where jt(s) is the spherical Bessel function of order 
We obtain 



y ' 1=0 



(48) 



where S g e is the observed galaxy transfer function [anal- 
ogous to the CMB radiation transfer function gi(k)\. 

$° g ,t = i l J dyf(y)S g M^)D(y)ji(ky). (49) 

Similarly, we can write Sg # as 

/df n - 6 
dy-±S g , SD (VD(y)j^fi(ky). (50) 



Now we use Eq. (|34p to express 5 g ^ in terms of the over- 
densities generated by infiaton and curvaton fields: 



S° g4 = i e J dyf{y)D{y){[l + b 1 {k)]Sf ) 

+ [l + b 2 {k)]8f ) ) n {ky) (51) 



and 

SI 



dy^-D(y)(sl 1) +5l 2) )^f l (ky). (52) 



Using these expressions together with equation (f36j) it 
is straightforward to calculate angular power spectrum 
Ci, which can be conveniently divided into three terms: 



(53) 



corresponding to the pure real-space galaxy density con- 
tribution (gg), the redshift space distortion term (vv), 
and the cross-correlation (gv). These components of the 
angular power spectrum are expressed in terms of the 3D 
linear matter power spectrum A 2 and x±, x 2 parameters, 
introduced in the previous section: 



4, / ^A2(|^°(fc)| 2 (l-x 2 ) + |^ 1 (fc)| 2 2 ; 2 ), 



C f = 8nJ ^At(M[wt(k)WZ(k)](l-x 2 ) 

+ U[Wl*(k)W[{k)]x 2 ), and 
CT = 47r| ^Al\W[(k)\ 2 , 
where the window functions are given by 
W?(k) = / bD(y)f(y)j e (ky)dy, 



(54) 



Wj{k) = I (b + Ab)D(y)f(y)j e (ky)dy, and 

2/ 2 4- 2/ — 1 



n» m »(r)D(y)f(y) 



(2£- l)(2^ + 3) 



£(£-1) 



(2£-l)(2£+l) 

(l + l)(l + 2) 
(2£ + l)(2£ + 3) 



h-i{ky) 



3t+2{ky) 



dy. 



.u(ky) 



(55) 



Here b is the standard Gaussian bias, while A6 is a cor- 
rection that applies to contributions from the curvaton 
field only: 



A6 = 3g(6-l)4 fc2r( ^ (r) 



VI. CONSTRAINTS 



(56) 



To constrain x\ and x 2 parameters using large scale 
structure we use the code developed and first imple- 
mented in Slosar et al. (25|. We included the same 
data: the 5-year WMAP bispectrum x\ — 51 ± 30 
(la) di|; and the SDSS data (spectroscopic and pho- 
tometric luminous red galaxies, the photometric quasar 
sample, and the integrated Sachs- Wolfe effect via cross- 
correlation). We further included ancillary data to con- 
strain the background cosmological model and break de- 
generacies with the non-Gaussianity parameters [x\, x 2 ): 
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the CMB power spectrum [64-68j and supernova data 

The Markov chain results are displayed in Fig.|3]where 
the probability density distribution is plotted on the 
(xi,X2) plane. Dark regions show regions with the high- 
est likelihood and the contours outlines 68%, 95% and 
99.7% confidence levels. As the role of the curvaton 
field in the primordial density perturbation decreases, i.e. 
x 2 — > (£ — >• oo), the upper limit on ii decreases. This 
is because at small x 2 , LSS becomes much more sensitive 
to Xx, as one may see from the x 2 in the denominator of 
Eq. ([SI)]): the \W}{k)\ 2 term in Eq. {5JJ scales as x?/x 2 . 
In particular, if a local-type CMB bispectrum is ever ro- 
bustly detected (xi ^ 0) then the non-detection of excess 
large-scale clustering in SDSS would immediately set a 
lower limit on xi- 

VII. DISCUSSION 

This paper has extended the analysis of non- 
Gaussianity constraints into a two field inflationary mod- 
els. In most previous studies of non-Gaussianity it was 
assumed that primordial density perturbations were gen- 
erated either by inflaton field, in which case they are 
perfectly Gaussian, or only by the second field (for ex- 
ample curvaton) which contains quadratic part and gen- 
erates non-Gaussian initial conditions. It is important, 
however, to realize the possibility of an intermediate case 
where part of the curvature perturbation is derived from 
quantum fluctuations of the inflaton field, while an ad- 
ditional part is associated with a second field and con- 
verted to an adiabatic perturbation upon its decay. This 
results in a peculiar type of non-Gaussian initial condi- 
tion (which we may call "local-stochastic" since the field 
4>2 entering in the nonlinear term is correlated with but 
not identical to the linear potential) that is both observ- 
able and distinguishable from the curvaton-only "local- 
deterministic" or /nl form. This type of non-Gaussianity 
has two parameters: a nonlinear coupling coefficient /nl, 
and the ratio £ of inflaton to curvaton contributions to 
the primordial density perturbation spectrum. We con- 
nect these parameters with parameters characterizing in- 
flationary fields in Eqs. ([9]) and (JTUJ) . 

Using the power spectrum and bispectrum constraints 
from SDSS and WMAP we are able to constrain these 
parameters. Adding two sets of constraints together al- 
lows us to break the degeneracy in the (/nl,£) param- 
eters that exists with the CMB bispectrum alone. If 



non-Gaussianity in the CMB is ever detected, and the 
bispectrum has the local configuration dependence, this 
will enable us to measure the relative contributions of the 
inflaton and curvaton. 

We have found that in contrast to the local- 
deterministic non-Gaussianity, whose main effect on 
large scale structure is a scale-dependent increase in the 
bias, the local-stochastic non-Gaussianity can introduce 
stochasticity between the matter and halo distributions. 
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FIG. 3: Constraints in the (xi,X2) plane, including both the 
CMB bispectrum and the galaxy power spectrum. 



It can also lead to relative stochasticity between haloes of 
different masses, since Eq. (j3"9")l depends on the Gaussian 
bias 6 g of the haloes (e.g. % — > 1 if fe g -> 0). The potential 
use of these effects to directly test the hypothesis of mul- 
tiple fields contributing to the primordial perturbations 
is left to future work. 
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